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It is well known that a closed universe with a minimally coupled massive scalar field always 
collapses to a singularity unless the initial conditions are extremely fine tuned. We show that the 
corrections to the equations of motion for the massive scalar field, given by loop quantum gravity in 
high curvature regime, always lead to a bounce independently of the initial conditions. In contrast 
to the previous works in loop quantum cosmology, we note that the singularity can be avoided even 
at the semi-classical level of effective dynamical equations with non-perturbative quantum gravity 
modifications, without using a discrete quantum evolution. 
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Various models in cosmology predict a dead end of 
our universe in a final state of collapse known as big 
crunch. It is well known that a big crunch may occur 
if the universe is matter dominated and starts collaps- 
ing under its own weight, leading to a cosmic doomsday 
scenario. The features of a classical cosmological evo- 
lution strongly depend on the sign of spatial curvature. 
In standard cosmology a flat or open universe can ei- 
ther experience only expansion or contraction without a 
transition from one regime to another, unless it is filled 
by some exotic matter with a negative energy density. 
Thus, if such a universe is expanding it never dies in a 
big crunch. On the contrary, a positive spatial curvature 
allows transitions from expansion to contraction and vice 
versa, moreover, in some models such transitions are in- 
escapable. With the current astronomical observations 
constraining VLq = 1.02 ±0.02 m, it is quite possible that 
our universe may be closed and heading towards a big 
crunch. Apart from standard cosmology, big crunch sce- 
nario is also present in models of string cosmology ^] and 
braneworlds 

In classical cosmology a k = 1 universe if filled with 
a matter source satisfying equation of state u) < —1/3 
(where uo = p/ p) could have a finite minimum size, since 
in this case the curvature term in Hubble equation dom- 
inates over the matter term for small a. Such a uni- 
verse will neither have a big bang nor a big crunch. On 
the other hand if w > —1/3 then the universe starts 
from a big bang and ends in a big crunch. The prob- 
lem of big crunch has been extensively studied in k = 1 
FRW model with a minimally coupled massive scalar field 
0, IE S Hi ■ Such a study captures the essential fea- 
tures of the dynamics before big crunch in general. Since, 
the effective lo for a scalar field cj) with some nonzero po- 
tential V{(t)) can vary in the range [—1, 1], hence a variety 
of possibilities exist. If V{(f)) is a smooth function pos- 
sessing a minimum (this type of potential is favored by 
a theory of reheating), then the scalar field at late epoch 
of cosmological expansion will mimic an ordinary mat- 
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ter with lj > 10] and the transition to a contracting 
phase proceeds. We should also notice that a wide family 
of potentials without a local minima invoked to explain 
issues of dark energy (see 0, 0, fo'" reviews) lead 
to eternal expansion even in the case of k = 1 , and such 
models do not face the big crunch problem. This is also 
true for some potentials with local minima as has been 
shown in |Q| & ^^1- In this work we will, however only 
consider potentials with a local mimina in a framework 
which leads classically to a big crunch. 

An important question which thus arises is whether 
a collapsing universe with a scalar field experiences a 
bounce or ends in a big crunch. Answer to this depends 
not only on potential V{4>), but also on a particular phase 
trajectory chosen. Bounce is indeed possible for poten- 
tials which are not very steep, but detailed investiga- 
tions have shown, that this in general requires a severe 
fine-tuning of initial conditions. Moreover, the revealed 
chaoticity of such cosmological dynamics indicates, 
that a bounce in the end of one cosmological cycle reflects 
nothing about a possibility to have a bounce during the 
next contraction stage of the universe. To ensure two 
consecutive bounces it is therefore necessary to further 
constrain the initial conditions. As a result only a zero 
measure set of initial conditions gives us a possibility to 
have infinite number of bounces, and then to escape a 
singularity for an infinite time interval. This leads to a 
conclusion that the future singularity problem can not be 
solved in classical cosmology with a massive scalar field. 

Occurrence of singularities in classical cosmology has 
always been thought of as a signature of a domain where 
general relativity ceases to be operational and must be re- 
placed by a quantum theory of gravity. In the present pa- 
per we try to understand dynamics leading to big crunch 
in light of one of the approaches to quantize general 
relativity known as loop quantum gravity. It is based 
on canonical quantization of general relativity and its 
salient features include background independence, non- 
perturbativeness and pred iction of discrete spectrum for 
geometrical operators 16]. This theory has recently been 
applied in cosmological scenarios and has yielded various 
novel results which include absence of big bang singu- 
larity [l7j, insights on initial conditions of the universe 
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[Tsl l , possibility of inflation with right e- foldings , 
suppression of power at large scales in CMB .2QJ etc. For 
related works in this direction see for example |2ll |. 

One of the key features of loop quantum gravity is that 
it predicts a modification in the behavior of geometrical 
density at short distances. Unlike in the conven- 
tional quantum cosmology where the geometrical density 
blows up when a — > 0, in loop quantum cosmology the 
eigenvalues of the corresponding operator become zero 
at a = 0. We should note that near a ~ the picture 
of a continuous spacetime breaks down and evolution is 
described by difference equations, however, there exists a 
domain where the continuous spacetime picture is valid 
and quantum gravity effects are manifest. The eigen- 
values of geometrical density operator in loop quantum 
cosmology are given as ,17], 

d,{a) = i?(aVa')a~' , a' - 7^1 j/3, (1) 

where a is the scale factor, £p\ is the Planck length, 7 
is the Barbero-Immirzi parameter and j is a half-integer 
greater than unity. Here the function D{a^ /a^) is derived 
from quantum theory and is given by, 

D{q) = {8/77fq^^^7[iq + 1)"/^ - \q - l|ii/4] 

- nq[{q + 1)^/4 - sgn {q - l)\q - 1]'/*]}^, (2) 

where q := a^/al. It should be noted that the scale 
where quantum effects become significant is defined by 
a* and though the Barbero-Immirzi parameter can be 
fixed to 7 = ln2/(-\/37r) by black hole thermodynamics, 
the parameter j is arbitrary and can be used to set an 
effective quantum gravity scale in the theory which may 
even be bigger than ^pi. 

The density operator in eq. has various peculiar fea- 
tures. In the quantum domain when a ^ a* it behaves 
as 

dj il2/7)^{a/a^y'^a-^ . (3) 

In transition to the classical regime when a increases, 
dj effectively becomes proportional to lower powers of a 
and eventually becomes dj ~ a^^ when a >> a*, which 
signals onset of a classical regime. 

The Hamiltonian of a massive scalar field (f> with a po- 
tential V{4i) in loop quantum cosmology can be written 
as 

n = d,^+a'V{^), (4) 

where is the canonically conjugate momentum to (f> 
and satisfies = d~^(j). For the case of k = 1 FRW 
model the effective Friedmann equation becomes 



together with the modified Klein-Gordon equation 

■ d Q 
6 — a H 4> — \tl dj — a dj V a, 
da 

= (^-3H+^^^~DV^ (6) 

which yield the Raychaudhuri equation with quantum 
gravity correction as 

d 2tt f I -.-i 1 -o d , \ Stt ,, 

a - -3M^ [d^ -^^ d-a'O^SMf/'^^^ 

= fl--^V^n0)- (7) 

Before we study some interesting phenomenon due to 
quantum gravity modifications in above dynamical equa- 
tions, it is important to emphasize the semi-classical na- 
ture of these equations. As stressed earlier, the picture 
of a continuous spacetime breaksdown near a ~ 0, where 
the description is entirely in terms of spin network states. 
Hence, the dynamical equations H5I6I7I) are not valid very 
close to a ~ 0. This is the regime where full quantum 
gravity is functional and the evolution can be described 
only in terms of quantum difference equations. The dy- 
namical equations are thus semi-classical in nature which 
are valid in the domain around a* , where we have a con- 
tinuous spacetime but with quantum gravity modifica- 
tions. It is also interesting to note that the full theory of 
loop quantum gravity does not have a notion of external 
time which is present at the semi-classical level, so the 
notion of trajectories which are solutions of the dynam- 
ical equations is invalid near a ~ 0. Thus if we wish to 
resolve the initial big bang singularity then one has to 
do a full quantum gravity treatment based on discrete 
form of the Wheeler-DeWitt equation However, as 
we would show, to avoid a big crunch singularity, the 
modifications of dynamical equations at the above semi- 
classical level are sufficient. 

Using eq. ^ it is easy to see that the ((> term in Klein- 
Gordon equation changes its sign as the universe evolves 
from a quantum regime to a classical regime or vice versa. 
This behavior is crucial in establishing various key results 
in loop quantum cosmology. For example, if the universe 
is expanding from a quantum regime the usual friction 
term in the Klein-Gordon equation is replaced by an anti- 
friction term which is responsible for an inflaton to climb 
the potential hill, even if it starts from its bottom, as 
has been shown for the flat model (19i. .20 ] and the close 
model 22]. It is important to note that the potential 
term in Klein-Gordon equation becomes sub-dominant in 
quantum domain, since the function Z) <C 1 for a <C a*. 

If we use the conventional Klein-Gordon equation with- 
out quantum correction in the case when universe is col- 
lapsing, then the 4> term in classical regime acts like an 
anti-friction term and the scalar field acquires very high 
values just before big crunch. However, in loop quantum 
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cosmology this situation will be prevented since when a 
collapsing universe reaches a size comparable or less than 
a*, the otherwise anti-friction term becomes a frictional 
term with a higher magnitude which stops the scalar field 
and leads to a bounce. In the quantum regime, the strong 
dominance of the kinetic term over the potential term 
means not only that the mechanism is robust to a change 
in V{(j)), but also that it dominates over gradient terms 
in scalar field. 

To demonstrate the effect of loop quantum gravity on 
big crunch let us first consider the case of a quadratic 
potential 

Vi^)^^mlcj,' . (8) 

In the classical regime the equation of a massive scalar 
field © describes the behavior of a harmonic oscilla- 
tor. A detailed analysis of cosmological solutions in this 
case can be found for example in jj]. It is important 
to note that in the classical regime when the (j) term is 
anti-frictional during collapse, almost all initial condi- 
tions lead to (possibly after a period of growing scalar 
field oscillations) a regime where the potential term in 
Friedmann equation (jsj is negligible in comparison to 
the kinetic term. In this domain, 

and a singularity is reached at some moment t = to, 
where to and c are constants of integration. Some tra- 
jectories can avoid this regime and experience bounce, 
giving rise to an interesting fractal structure of the phase 
space (see 0, for details) , but their measure is very 
small. Such a possibility is approximately inversely pro- 
portional to amax/ip\, where amax is the scale factor at 
the point of maximal expansion U . 

Solving cqs.|(Sl-|7|) along with eq.(^ numerically we can 
understand the change in dynamics due to quantum grav- 
ity effects and as expected, the change from the classical 
equations is very drastic. The evolution of scale factor 
for the quadratic potential is shown in Fig. We con- 
sider initial conditions such that the universe encounters 
a big crunch classically (as shown by a dashed curve). 
The evolution of scalar field with respect to scale factor 
is shown in Fig. [2]. As can be seen near the big crunch 
the field (p takes very large values. It is clear from both 
of these figures that this situation is averted in quantum 
gravity. 

When the size of the universe becomes comparable to 
a, and less, the quantum corrections become important 
in dynamical equations. As mentioned earlier, in the 
collapsing case the anti-friction term which was —iH(j) 
(with negative H) in classical regime, changes sign, be- 
comes frictional and even goes as 12H(j) for a ^ a,. In 
this process, quantum gravity essentially applies breaks 
to the motion of through a large friction term in Klein- 
Gordon equation and the scalar field almost completely 
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FIG. 1: Evolution of scale factor for a closed universe in classi- 
cal cosmology is shown by the dashed curve and in loop quan- 
tum cosmology is shown by the solid curve for the quadratic 
potential with j — 100, = O.lMpi and initial conditions 
ai = 10 a,, Hi =0,(t>i — and (j>i determined from eq.Q. A 
collapsing closed universe in classical cosmology encounters 
a big crunch whereas in loop quantum cosmology it bounces 
into an expanding phase when a — 0.38 a, and avoids big 
crunch. 

"freezes" at some high value (pf. At the instant when the 
motion of cj) ceases and (f) — 0, all quantum gravity signa- 
tures are excluded from the model and equation of state 
becomes —1. Such equation of state in the case of posi- 
tive spatial curvature ultimately leads to a bounce. So, 
the contracting universe bounces and turns to expansion. 

After the bounce when the universe is still smaller than 
a*, the (f) term in eq.© acts like an anti- friction term 
(since, the Hubble parameter H changes sign to positive 
in expanding phase). The scale factor in this phase grows 
very rapidly and eventually as a 3> a* , the quantum cor- 
rections become negligible and the (f) term becomes fric- 
tional. The universe then enters a standard infiationary 
stage. An example for this behavior can be seen in vari- 
ation of scale factor in Fig. After checking a large 
range of initial conditions and parameters j and we 
found that a bounce always occurs and the behavior of 
dynamics is qualitatively the same as shown in Figs. ^ 

With the occurrence of bounce established for the 
quadratic potential in a generic way we now examine the 
case for a steep potential, 

V{cP) = 2 (cosh(0VMp2 ) - 1) . (10) 

The dynamics of a scalar field in classical cosmology is 
rather different for potentials less and more steep than 
the exponential one. A steep potential never becomes 
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FIG. 2: Plot of (j) in Planck units and a/a, in the domain 
where the transition from classical to quantum and then again 
to classical occurs in loop quantum cosmology for 0^ poten- 
tial. The dashed curve is the one for classical cosmology and 
the field approaches infinity as universe encounters big crunch. 
The solid curve is for loop quantum cosmology evolution. In 
this case it is seen that the field freezes near the bounce. Ini- 
tial conditions, j and are same as in Fig. 0. 
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FIG. 3: Scale factor for the case of steep potential llOt at 
late times (in Planck units) of collapse. Initial conditions are 
ai = 100a*,iifi = and (fti = Q (which also determine cj>i) 
with j = 100. The evolution in loop quantum cosmology is 
depicted by solid curve which avoid big crunch when a = 
0.082a* whereas classical cosmology shown by dashed curve 
encounters a singularity. 



dynamically unimportant at a contraction stage, and a 



scalar field oscillates infinitely while falling into a singu- 
larity |2^. The bounce for steep potentials as classically 
understood is impossible and a thus a singularity is in- 
evitable 0. 

These features, however, do not change the dynamics 
in quantum regime. The quantum gravity corrections 
are significant mainly due to their influence on terms 
in dynamical equations and thus in the quantum regime 
they become dominant over the potential terms. The nu- 
merical results for the potential UlUfl are shown in Figs. 
PI & 2] which clearly show a bounce. The behavior of 
the universe near bounce is very similar to one shown in 
Figs. ^&[2]. In Fig. 0] we have shown how in clas- 
sical regime the field undergoes oscillations and then 
in quantum regime it first freezes and bounces back. As 
for the quadratic potential we found that the bounce for 
steep potential (|10(l always occurs for arbitrary choice of 
initial conditions and j. 
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FIG. 4: Dynamics of (j> (in Planck units) and a/a, for po- 
tential (IIOI I with initial conditions same as in Fig. j^. In 
the classical regime the field undergoes oscillations. In loop 
quantum cosmology the field freezes near bounce (shown by 
solid curve) whereas in classical cosmology (shown by dashed 
curve) it takes large values near big crunch. 

Our results are in a sharp contrast to earlier results 
in classical cosmology, where a occurrence of a bounce 
requires a very special initial condition to be chosen. As 
the quantum bounce is independent of initial conditions, 
it eventually occurs also in the second, third and con- 
secutive cycles of contractions giving rise to some kind 
of aperiodic cyclic universe. It is interesting that the 
bounce via loop quantum gravity is rather similar to an 
artificial bounce put by hand into the classical picture 
in j2j| (d — > —a and (p ^ <j) a,t some scale factor of the 
order of the Planck length). This suggests that analysis 
of effects of "hysteresis" of the scalar field evolution and 
time asymmetry which have been done in |24| are also 
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applicable to our model. 

We should stress that for both considered potentials, 
quadratic (jHJ and steep ifinjl . we found that the bounce 
occurs when size of the universe is not very small com- 
pared to . For the quadratic potential a closed universe 
bounces when a ^ 0.2 — 0.5 a* and for a steep potential 
it bounces when a ^ 0.05 — 0.3a,. We also found that 
the role of j on the nature of dynamics near big crunch 
is very weak. It is important to note that the big crunch 
singularity has been avoided purely at the semi-classical 
level of effective dynamical equations with quantum grav- 
ity modifications. This is in distinction to the earlier 
work in loop quantum cosmology where the singu- 
larity avoidance was accomplished using discrete quan- 
tum evolution near a = 0, in which case one can not use 
the picture of a continuous spacetime and the dynamical 
equations are invalid. 

Thus, we establish that for a closed universe with a 
minimally coupled massive scalar field a bounce will oc- 



cur for any choice of initial conditions as soon as uni- 
verse becomes of the size that its dynamics is governed 
by equations with loop quantum gravity modifications. 
Hence a big crunch or a cosmic doomsday will always be 
avoided. In future, it will be important to study dynam- 
ics of bounce in anisotropic case and also the role of mat- 
ter and radiation near bounce. Apart from loop quantum 
cosmology, bouncing universe scenarios have been in dis- 
cussion in various string models (for example see |25l|'). It 
will be interesting to investigate their common features 
with loop quantum effects for various values of geomet- 
rical curvature. 
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